Enhancement of electron-phonon energy transfer and ultrasound attenuation in 
impure conductors due to Mandelstam-Leontovich relaxational mechanism 
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We derive a general relation between the ultrasound attenuation rate r~ h (w) due to interaction 
with conduction electrons and the energy flow J(T) between overheated electrons and phonons at 
bath temperature (cooling power). We demonstrate an existence of several mechanisms leading to 
strong enhancement of T~ h (w) and J{T) at low frequencies/temperatures in disordered conductors; 
all of them are related with the Mandelstam-Leontovich relaxation mechanism. The ultrasound 
attenuation and the electron-phonon energy transfer are found to be sensitive to magnetic field and 
can be enhanced by the factor 10 2 — 10 3 at B ~ 10 Tesla. They can also be increased strongly in 
semiconductor heterostructures with an external gate. 



The theory of electron-phonon interactions in disor- 
dered conductors has a long history. One of its key 
points, known as "Pippard ineffectiveness condition" 
(denoted as PIC below) [1, 2], is that for phonons with 
small momenta q the rate of inelastic electron-phonon 
processes is by factor ~ ql <C 1 weaker than in the ballis- 
tic limit (/ is the elastic scattering length of electrons) . In 
the framework of this theory the ultrasonic attenuation 
rate was calculated [1] in the ql -C 1 limit: 
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where v is electron density of states [4], pf = itivf is 
Fermi momentum, n e and p rn are the density of elec- 
trons and the mass density of material, D = vpl/d e 
is the diffusion coefficient, q is the phonon momentum, 
d e is the dimensionality of electron motion. The sub- 
script a corresponds to the choice of either transverse 
(tr) or longitudinal (1) phonons; correspondingly, numer- 
ical coefficients c a arc defined as c tr = 1/(2 + d e ) and 
q = 2(1 — d ( T 1 )/(2 + d e ). In the clean limit ql » 1 cou- 
pling to transverse phonons is ineffective, whereas longti- 
tudinal phonon attenuation rate scales as oc q. 

Below we will see that phonon attenuation rate 
is directly related to the energy flux (cooling power) 
J{T e ,T p h) = J{T e ) — J(T p h) transmitted from electrons 
with temperature T e > T p h to phonons. As the tempera- 
ture decreases below To = Hv s /l, wave vectors q = T/hv s 
of thermal phonons drop below I -1 and crossover takes 
place from the classical result for the e-ph cooling power 
(assuming bulk phonons) J(T) oc T 5 , to the faster de- 
crease with temperature J(T) oc T 6 /T . The latter be- 
havior of the cooling power predicted in Ref. [2, 3] cor- 
responds to the result of Ref.[l], Eq.(l) for the phonon 
decay rate, and was shown to be obeyed very accurately 
by the data on disordered metal films of Hf and Ti [7] . 

It is known that PIC does not hold if either a) im- 
purities do not follow the lattice deformations [9], or b) 
piezoelectric interaction is present [8]; enhancement of 



the cooling power in semiconductors due to the "electron- 
phonon-impurity" interference has been also discussed 
[10]. Data on the cooling power in some strongly dis- 
ordered conductors [11-13] where none of the conditions 
a) or b) are expected to hold, indicate to the presence of 
a different enhancement mechanism. 

In the present Letter we demonstrate an existence of 
quite general mechanism for the enhancement of both 
cooling power J(T) and ultrasound attenuation rate 
T ph (f° r longitudinal phonons) at low temperatures. 
This mechanism is similar to the Mandelstam-Leontovich 
(ML) relaxation [14] and is effective if lattice motion is 
able to induce significant oscillations of local densities of 
certain globally conserved physical quantities. The devi- 
ations of these local densities from their equilibrium val- 
ues are enhanced by slow diffusive character of electron 
motion (characterized by both small frequency and small 
momentum) aimed to restore equilibrium. This leads to 
a significant retardation in the response and thus to the 
entropy production and dissipation. In contrast to PIC 
which suppresses the relaxation rate at strong disor- 
der and small carrier concentration, the ML mechanism 
is efficient at these conditions. In particular, we found 
a broad range of frequencies/temperatures where r^ 1 is 
^-independent, whereas J(T) oc T 4 . 

We will discuss below two specific realizations of such a 
mechanism: 1) the presence of several types of quasipar- 
ticlcs with different spectra which can be described by the 
quasi-spin components, and 2) relatively weak Coulomb 
interaction between electrons, when the local electroneu- 
trality condition is not strictly obeyed. In the first case 
phonons create non-equilibrium distribution of quasi-spin 
density, while in the second case it is the charge den- 
sity. We will show that the physics behind calculations 
in Ref. [10] is precisely the charge density relaxation at 
incomplete screening rather than the peculiarity of the 
momentum-averaged e-ph vertex. 

We will use the co-moving frame of reference (CFR) 
rigidly bound to the lattice. Then electron-phonon inter- 
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action is encoded mainly in non-uniformity of the coordi- 
nate transformation. For a single branch of electrons, the 
density dependence of the Fermi-energy leads to[l, 15]: 

H e -ph = y]p a (vp'Vpu a ) q ■ -Ipp^p+q (2) 
p,q 

where p and v denote electron momentum and velocity, 
respectively and u is the lattice displacement. The tensor 
structure of Eq.(2) is crucial for local processes only, such 
as those giving rise to the result given in Eq.(l). To study 
the enhancement of dissipation due to diffusion processes, 
it is sufficient to average the e-ph vertex over the Fermi 
surface. For a metal with spherical Fermi-surface one 
finds PaivpV pu a ) q = Tdivu, where Y = ppVF/d e . In 
general T may contain other contributions. In particu- 
lar, for semiconductors T is known [16] to be much larger 
than Ep . Under the condition of strict electroncutrality, 
the scalar vertex T is screened out completely and the 
classical result Eq.(l) is valid. This is not the case, how- 
ever, when N different types of quasiparticles are present. 
Then the interaction can be written as 

N 

H e - ph = J2 Ij divu (3) 

i=i 

The Coulomb interaction is able to screen out only the 
single mode corresponding to the total density n = 
^2j(iJip)j, whereas N — 1 asymmetric modes survive 
screening [15]. In particular, for the N = 2 case cor- 
responding to the Zecman splitting of spin-up and spin- 
down sub-bands, the mode corresponding to spin density 
fluctuations is unscreened. As we will see below, coupling 
to the spin density mode generates diffusion motion of 
electrons which leads to extra dissipation seen in the in- 
crease of both the ultrasound attenuation and of the heat 
transfer from hot electrons to the phonon bath. The same 
effect occurs even in the case of a single branch of quasi- 
particles if Coulomb screening is incomplete. The lat- 
ter is natural for 2D (or quasi- ID [10]) electron systems 
where Coulomb interaction mainly propagates through 
the surrounding dielectric media. 

We start by writing down the quantum kinetic equa- 
tion in order to establish a general relationship between 
the rates of electron cooling and ultrasonic attenuation. 
In a typical experiment on electron cooling the elec- 
tron system possesses its own temperature T e i decoupled 
from the lattice with the temperature T p h- Assumption 
of a well defined electron temperature is reasonable if 
electron-electron inelastic rate is much larger than 
electron-phonon one, T~_ ph , which is the case at low tem- 
peratures. Under that assumption of local equilibrium, 
electron-electron interaction drops out of the theory. Our 
results for r~_ h are presented in [17]. 

The key point to establish a relationship between elec- 
tron cooling and phonon relaxation, or ultrasound atten- 
uation, is to consider the kinetic equation for phonons 



instead of the electron kinetic equation and also to as- 
sume conservation of the total energy of the electron and 
the phonon systems. Thus we consider phonon heating 
rather then electron cooling, which makes it possible to 
express heat flow in terms of equilibrium electronic quan- 
tities (see [17], Sec. II): 

d t B ph (uj,T ph ) = [B ph (ui,T e i) - Bphiuj^Tph)}-^ (u,T el ) . 

(4) 

where B ph (uj,T) = B ph {u/T) = i(coth(w/2T) - 1) is 
the phonon distribution function. We employed here 
the relation t - , 1 = —=— [imS^fu.q)! between the 
phonon life-time T p ^ and the retarded phonon self-energy 
£ fl [17]. Now the energy flow J = can be written 
as follows: 

oo 

t f j < \ Bphiu/Tet) - B ph (uj/T ph ) 

J = / duu^ uP F — (5) 

J T p h(0J,T e l) 


where u p h(uj) = u 2 /(2tt 2 v^) is the phonon density of 
states for 3D phonons with the sound velocity v s . 

One can see from the relation Eq.(5) that the cool- 
ing power is determined by the phonon relaxation rate 
T p h(ijJ , T e {) which also determines the ultrasound atten- 
uation. If we do not take into account any correlations 
induced by electron-electron interactions (charge screen- 
ing is considered in the RPA approximation only), then 
T p h does not depend explicitly on electron temperature 
T e i allowing to separate clearly the incoming and outgo- 
ing energy flows: J — J(T e {) — J{T p h). In particular, 
for the power-law w-dependence of T p ^ 1 (w) oc w' 3 the in- 
coming energy flow for the 3D phonons is proportional 
to J{T el ) cx T e 4 ; +/3 , for details sec [17], Sec.V. 

Below we present a simple derivation of the diffusion- 
enhanced contribution to the phonon relaxation rate T~ h 
in terms of macroscopic equations for the current and 
density of electrons; alternative diagrammatic derivation 
is presented in [17], Sec. III. In the CFR the continuity 
and diffusion equations for each i-th species of quasipar- 
ticles read: 

d t n (i) +divj w = 0, 

where (i) stands for the quasiparticle branch number, riW 
is the electron density, is the particle number current, 
Ki = i>iDi is the mobility, Di is the diffusion coefficient for 
the i-th branch, = -VC/ (i) and U {€) is the potential 
energy. In the simplest derivation we assume no inter- 
branch scattering and thus the continuity equations in 
Eq.(6) imply that each of the partial electron densities 
n W ar e conserved separately. Generalization to the case 
where there is mixing between the branches will be done 
at the end of the paper. The potential energy 
Uc + in Eq.(6) consists of the Coulomb part Uc and 
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the phonon part $W = rW divu: 



V (r-r')^^(r') + T« cli^ 

3 



(7) 



where Vo(r) is the bare Coulomb potential acting be- 
tween conduction electrons; below we use its Fourier- 
transform Vo(q). Note that Vo(q) does not include screen- 
ing by conduction electrons in the sample. 

First, we discuss the case of several quasiparticle 
branches and perfect electroneutrality. The latter for- 
mally corresponds to the limit Vo(q) —> oo. For the den- 
sity modulation ny)(w,q) induced by the phonon with 
frequency u> and momentum q one finds from Eqs.(6),(7): 



(w,g)=n i (w )9 )(* i (w J «)-# c (w, g )) (8) 



where q) = divu , $ c = £\ Ej n j rep- 

resents dynamical screening of Coulomb interaction and 
Ilj = ViDiq 2 1 '(— icj + Diq 2 ) is the polarization function. 

The solution (8) obeys charge-neutrality: n to t = 
J2 n ^ — 0- However, there arc (N — 1) independent 

i 

density modes which are excited by phonon modula- 
tion. To simplify the expressions we consider below the 
case N = 2. 

The contribution of ML mechanism to the phonon de- 
cay rate may be expressed as t~£ ml = where Q t 
and E w are the dissipation power and the acoustic wave 
energy in a unit volume, respectively: 
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Here u m is an amplitude of ionic displacement and u = 
(q/q)u m exp[— iujt + iq ■ r]. Simple calculation leads to 
the following expression for the "diffusive" contribution 
to the decay rate of acoustic phonon: 
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where v s = uj/q is the sound velocity, while = (v 1 1 + 

are the effective density 



andD* = (D^+D 2 



i\— l 



of states and diffusion coefficient, respectively. 

In particular, application of parallel magnetic field H 
to a two-dimensional electron gas in a semiconductor 
hetero-structure leads to the splitting of spin sub-bands 
with r-f — Tj. = uH ([17]), where fi = gag is the elec- 
tron magnetic moment. Then, according to Eq.(10), the 
phonon relaxation rate acquires an 7J-dependcnt contri- 
bution that may become dominant at sufficiently strong 
field and low phonon frequencies. Adding the classical 
result (1) and the magnetic-ficld-controlled contribution, 
Eq.(10), one finds ([17], Sec. III.B, VII) for the full 
phonon decay rate = [t^]" 1 ■ Fh{<1i h), where for q 
parallel to 2D gas: 



T H (q,h) = 1 



{Dqf 



(11) 



Here is given by Eq.(l) for longitudinal phonons 
and h = /j,H/2ef (we assume here h to be rela- 
tively small compared to unity). The enhancement fac- 
tor Th can become very large for strong spin polariza- 
tion, h ~ 1. In particular, for low phonon momen- 
tum, ql < v s /vf, the factor Tu is of the order of in- 
verse adiabatic parameter (vf/v s ) 2 ~ 10 5 . The same 
situation may arise in multiple- valley semiconductors (if 
inter-valley scattering by disorder is weak for some rea- 
son), where splitting is present without magnetic field 
and usually Ti — T 2 ~ r^IV Another relevant ex- 
ample is provided by ferromagnetic metals with strong 
intrinsic band-splitting due to the exchange field. In 
the case of Fe: uH* m 1.8 eV, ef = 11.1 eV, vf = 
1.98 x 10 8 cm/s, v s sa 6 x 10 5 cm/s resulting in the maxi- 
mum enhancement of the phonon relaxation time as large 
as T H ~ {uH*/e F ) 2 {v F /v s ) 2 ~ 10 4 . 

In the case of a single quasipaticle branch, diffusion- 
enhanced dissipation may appear due to violation of 
the charge neutrality condition at a large screen- 
ing length. To calculate it we employ Eqs. 
(6,7) and obtain for the density fluctuation n = 

vDq 2 [-iuj + Dq 2 + 2k q 2 V (q)] _1 q), which leads 
to the following contribution to the relaxation rate 



T P h]c(<l) 



1 



vDq A 



p m u 2 + (Dq 2 + 2Kq 2 V a {q)f 



(12) 



Comparison of Eq.(12) with the result of classical theory 
(1) gives the enhancement factor J-c{q) in the form 
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{v s /v F ) 2 + de A (q 2 l 2 )(l + 2vV Q (q)) 2 



(13) 



For the 2D gas with Coulomb interaction and the con- 
stant dielectric permittivity s of surrounding media we 
have Vo(q) = V2D(q) = 2-Ke 2 /eq. In the relevant range 
of q parallel to the 2D gas the Tc (q) factor reduces to a 
constant (similar result was obtained in [10]): 
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cj 1 fehv F \ 2 f r y 
9n \ e 2 J \PfVfJ 



(14) 



where g\j is dimcnsionless conductance per square in e /h 
units. Eq.(14) is valid both for 2D (c^ 1 = 4) and for 
quasi-2D (c^ 1 = 15/4) systems as long as the phonon 
wavelength A is much larger than the width of the quasi- 
2D system. Note the e 2 dependence of the ^-factor on 
dielectric susceptibility e. 

Even larger enhancement arises in the presence of a 
gate that additionally screens Coulomb interaction and 
allows the density to fluctuate stronger. For this geome- 
try and q parallel to 2D gas Vo(q) — > V g (q) = V^oiq) (1 — 
e~ 2gb ), b being the distance between the 2d electron gas 
and the gate. For phonons with the wavelengths l/q>b, 
the effective potential V g (q) « V(q) ■ 2qb w const and the 
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presence of adiabatic parameter in the denominator of 
(13) does become important at low enough temperatures: 
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16(T/ PF v F ) 2 



A(v s /v F ) 2 + (q 2 l 2 )(8niye 2 b/e 



(15) 



where Coulomb interaction is still assumed to be rela- 
tively strong: 2-Kve 2 b/e ^s> 1. 

Note that Eq.(15) allows to determine the background 
dielectric constant e in a conducting sample(previous ex- 
ample of screening due to an external gate can be con- 
sidered specific example with dielectric function 
e(q) = e/(l — e~ 2bq )). This part of dielectric function is 
formed by relatively high-energy processes and is usually 
added as a phenomenological constant in the theory of 
electromagnetic response of impure metals and weak An- 
derson insulators. In band insulators and semiconductors 
e can be measured spectroscopically by studying spectra 
of F-ccntcrs and Wannier-Mott excitons. In impure met- 
als spectroscopic methods fail due to strong absorbtion 
by conducting electrons while the direct measurement 
of static polarizability gives mainly the contribution due 
to low-energy transitions in the conduction band. One 
can see from Eq.(15) that the crossover frequency u cross 
where t~^uj~ 2 start to decrease as a function of u, is pro- 
portional to e: ui cross T = (v s /v F ) 2 (47ri/e 2 &) -1 e, where 
t = l/vp is the elastic scattering time. Thus measur- 
ing this crossover frequency one can estimate the back- 
ground dielectric constant e which controls the strength 
of Coulomb interaction and which is important for the 
theory of strongly disordered superconductors [18]. 

Finally, we collect results of both the ML enhancement 
due to the charge density and the spin density fluctu- 
ations, taking also into account mixing of spin projec- 
tions by the spin-orbit interaction characterized by the 
parameter r^ 1 . We also consider the dependence of re- 
laxation rate on the direction of phonon propagation rel- 
ative to 2D gas[17]. Both effects lead to the replacement 
Dq 2 /(-iuj + Dq 2 ) => Dq?,/(-iu + Dq 2 + r" 1 ). It results 
in the total enhancement factor of the form: 



T = Tc 2 



%q\v F h 2 



(qv s ) 2 + (Dq 2 + 1/t so ) 2 



(16) 



For 2D electrons and 3D phonons |q||| = q smQ is the 
phonon momentum component parallel to the 2D sys- 
tem which appears in all the terms originating from elec- 
tron diffusion. In this case J~c 2 d i s independent of q, 
and T has a maximum as a function of uj. The spin 
fluctuation effect given by the second term vanishes at 
small lu because of the mixing of branches caused by 
spin-orbit interaction. It also decreases at large u be- 
cause the dissipation power increases slower with ui than 
does the acoustic wave energy. At large enough Zee- 
man splitting h when the effect of spin fluctuations in its 
maximum is large, there is a wide frequency region (the 
falling part of the curve J- oc u;~ 2 in Fig.l) where 



is almost frequency independent. In this region the cool- 
ing/heating rate J(T) cx T 4 InT for the quasi-2D case. 
This temperature dependence is almost the same as in 
the case of impurities which are not fully involved in the 
lattice motion [9] . The extra logarithmic factor arises be- 
cause of the angular averaging of 1/T p h(0) dominated by 
the small values of 8. To illustrate this behavior we con- 
sider a thin film of semiconductor InSb (g- factor g w 50). 
At strong (and parallel to the 2D plane) magnetic fields 
gjJLBH Aso classification in terms of the spin sub- 
bands is still valid approximately, in spite of the Rashba 
spin-orbit coupling Aso- The analysis presented in [17], 
Sec. IV, VI, leads to Eq.(16) and is summarized in Fig.l. 
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(Color online) The total enhancement factor T = 
7r/2 of ultrasound attenuation in the 2D semi- 
conductor InSb. The parameters taken are n = 10 11 cm -2 , 
PfI = 50, Aso = O.lmeV and magnetic fields are 3T (blue), 
5T (red) and 7T (green). Dashed curves represent the result 
in the absence of SO relaxation, Aso = 0. 

In conclusion, we presented an intrinsic relation be- 
tween electron-phonon cooling and ultrasound attenua- 
tion due to conduction electrons and demonstrated an 
existence of a general relaxation mechanism that leads 
to enhancement of both the e-ph cooling power and the 
phonon decay rate. In particular, this relaxation mecha- 
nism is important for semiconductor heterostructures due 
to weakness of Coulomb interaction between conduction 
electrons; measurements of ultrasound attenuation may 
be used to determine the background dielectric constant 
s. Another version of the same mechanism is operative 
in the presence of weakly mixed electronic sub-bands. In 
particular, it may lead to a strong dependence of the 
cooling power on the external magnetic field or on the 
intrinsic exchange field in ferromagnetic metals. 
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Supplementary Material 



Keldysh component (...) may be parametrized as 



D K = D R oF-Fo D A 



(S19) 



I. HAMILTONIAN 

We consider an electron system with several spectral 
branches (which are labeled by i £ (1, iV) ). Interaction 
between electrons is considered in the direct (density- 
density) channel only, and will be treated within the ran- 
dom phase approximation (RPA). Electrons also scatter 
off local impurities, and we assume this scattering to be 
identical for all spectral branches, Ui mp s{v) = Ui mp (r). 
While treating electron-phonon inelastic processes, it is 
convenient to work in the co-moving frame of reference 
(CFR), where the Hamiltonian acquires the following 
form [1-3] in the momentum representation: 

H = H e0 + i? e -e + He-ph, (S17) 



H e -e = r X] V r (q)n q n_ q , 



H e -ph = y^r.i(divu' 



q /fc — q,zj 



where Ui mp (r) is disorder potential, Vo(q) is a bare 
Coulomb interaction attenuated by the background di- 
electric constant e which could be additionally screened 



by a metallic gate; 



Ep VWq.iVW and n c 



are the i-th branch partial- and total densities, 



respectively. I\ and (divu) q = iqu are the e-ph interac- 
tion constant averaged over the Fermi surface and the di- 
vergence of ionic displacement field u. The Hamiltonian 
(SI 7) contains no inter-branch scattering. We assume 
these processes to be much weaker than the intra-branch 
scattering and discuss their role later on in Sec. Ill C. 



II. DERIVATION OF THE PHONON KINETIC 
EQUATION 

Here we derive the phonon quantum kinetic equation 
via the Keldysh diagrammatic technique (see Ref.4 for 
details). For simplicity, we consider only the longitudinal 
phonons; the transverse phonons may be analyzed in a 
similar way. The Green's functions are matrices in the 
Keldysh space: 



D 



D- 



D A 




(S18) 







S A 



where Do and D are the bare and exact phonon Green's 
functions, respectively and S is the self-energy part. The 



where F is related to the phonon distribution function; in 
equilibrium F = coth(aj/2T). The sign o means the con- 
volution in the time domain. A bare retarded(advanced) 
phonon Green's function is 



R(A) 



1 



1 



9 9 ' 

v s q 



(S20) 



where p m is the mass density. To derive the kinetic 
equation one employs the Dyson equation for the matrix 
Green functions which together with Eq.(S19)results in: 

(5 - 1 -s)5 = i, (S2i) 

=>■ [Dq\F] = ^-S H of-f oE A (S22) 

According to Eq.(S20) Dq 1 contains the second deriva- 
tive with respect to the corresponding time argument of 
F = F(ti,t2). Then the commutator in Eq.(S22) reduces 
to the first derivative dt with respect to the "slow" com- 
bination t = t\ + 1% multiplied by the Fourier transform 
— iuj of the derivative with respect to the "fast" combi- 
nation t\ —ti. Averaging over the whole volume of the 
sample, one finds 



2ip m uid t F = Yi K — Yi R o F - f o E' 



(S23) 



In this equation we assume F{t\,t2) — > F(ui;t). We also 
omit the terms describing the external source that pumps 
energy into electronic system. 

A common action of the source term and the e-ph 
energy relaxation leads to a stationary energy distri- 
bution of both electrons and phonons. If the e-ph en- 
ergy relaxation is much slower than the electron-electron 
one, a quasi-equilibrium situation with two temperatures 
T e i, T p h is realized [5]. In this case the phonon distribu- 
tion function is F(uj,T p h) = coth(oj/2T p / t ) while Y, K is a 
quasi-equilibrium quantity. If in addition a weak e-ph in- 
teraction is assumed, electron-phonon interaction enters 
only trivially as the square of the e-ph matrix element 
in the Keldysh self-energy part S . In this approxi- 
mation Y* K depends only on the electron temperature, 
E* = J F(u;,T ei )(£*-£ A ), 



2i Pm ud t F{u,T ph {t)) = [F-(Z R 
-[F-(£«-£ A )] (u,q,T ph ). 



S A )] (oJ,q,T el ) 
(S24) 



Since the phonon decay rate is relatively low, r~ h <C u>, 
phonons are well-defined quasiparticles and the quasipar- 
ticle distribution function is sharply peaked around the 
phonon "mass-shell" u = sq. Thus, the quantities enter- 
ing in the R.H.S. of Eq.(S24) should be taken at oj = sq, 
see Ref. [4]. For convenience we introduce a standard 
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phonon distribution function B(oS) = (1/2) (F(uj) — 1) 
and define a decay rate T~ h : 



T P h(u,T e i) 



-ImS fl ( W ,q,T e/ )U = 



sq • 



(S25) 



d t B{u,T ph {t)) = B{^T el ) - B{u> T ph ) _ 

Tph{W,I e l) 

In order to obtain an electron-phonon heat flow, we have 
to multiply Eq.(S26) by both phonon density of states 
v p h{w) = w 2 /2ir 2 Vs and by energy, and integrate over u: 



Finally, let us mention that in 2D electron systems 
the phonon decay rate T p h(<l) depends on the angle 9 
between the phonon momentum q and the direction nor- 
mal to the plane of the 2DEG. Thus, in the equations 
like (S27,S28,S29) an angle-averaged decay rate should 
be used: 



')e H = \ f( shi9dB ) V>> 0) (S30) 



III. DIAGRAMMATIC DERIVATION 



J = I duv p h,{u)udtB(u),T v h(t)) 



(S27) 



dwv ph (uj) 



T p h(uJ, T e l) 



[B(u,T el ) - B(u,T ph )] 



The incoming and outgoing energy flows may be defined 
as follows: 

J {Tei , T ph ) = J+ (Tei ) - J- (T ph ,T el ), (S28) 

OO 

J+{T E i) = I dwv ph (u) 1 —= B(uj,T e i), 

J T p h{UJ,lel) 


oo 

J-(Tph,T ei ) = / du)v ph (u) B(uj,T ph ), 

J T p h{OJ,lel) 


In general, J_ depends on both T e i and T ph due to 
the effect of electron-electron interactions on the decay 
rate , see Eq.(S25). However such an effect is ab- 
sent within the RPA approximation for charge screening 
which we use here. Thus our result for the phonon life- 
time does not depend on temperature explicitly, t p h = 
T p h(y>). In this case the expressions for heat flows are 
simplified and J_ becomes a function of the phonon tem- 
perature T p h only. 

The electronic temperature relaxation rate t^ 1 can be 
obtained from the expression for the heat flow 3e- P h- 
Each phonon branch contributes as 



1 d^e—ph^a 

a ' 



dT PJ 



(S29) 



uv ph {u) 

2C e T J sinh 2 (w/2T) 
o 



where C e oc vT is the electronic specific heat. The full 
rate is then t^j = t^] + (d ph - l)r^ tr , where (d ph - 1) is 
the number of transverse phonon polarizations. Electrons 
are characterized by the quasi-equilibrium distribution 
function if their intrinsic inelastic scattering time r ee is 
much shorter than the temperature relaxation time te 
due to interaction with phonons. 





^ + v 



FIG. S2: Bare and screened e-ph vertices. Black, curvy 
red and dashed blue zig-zag lines are electron, phonon and 
Coulomb interaction propagators respectively. Bold zig-zag 
lines stand for statically screened Coulomb interaction (with 
screening by empty electron bubbles only). 



We present here the derivation of some of our results 
in the standard diagrammatic form, as it was done in the 
most papers on this subject [1-3]. 

The bare electron-phonon vertex corresponding to the 
Hamiltonian (SI 7) is of tensor structure in the space of 
electron species: 



V 



(S31) 



r 



N 



The diagonal structure of (S31) corresponds to our as- 
sumption about the absence of inter-branch mixing. Now 
one should take into account static Coulomb screening, 
which generates scalar counter-term 



■ 1 



(S32) 



The full screened vertex is then a sum T f u u = Tb are +Tc- 
The structure of these vertices is presented in Fig.S2. The 
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deformation potentials I\ averaged over FS are usually 
approximated as[l] 



r, = [(r 6s (p))| FS +p F v F /d e ] 



(S33) 



where Tb s represents the lattice-induced deformation po- 
tential under the lattice strain and p F v F /d e represents 
the averaged electron liquid stress tensor. We will con- 
sider the simplest model where the lattice contribution is 
uniform in momentum space Tb s (p) = Tt, s and thus is re- 
duced to the shift of electron band (see however Rcf. [6]). 




FIG. S3: a) Dynamically screened Coulomb interaction, 
where diffusion is taken into account. Namely, an impurity 
ladder is summed up b) The impurity ladder. Here black 
dashed line represents impurity correlator (£7(r)[/(r')). c) Di- 
agrams contributing to phonon self energy. This figure encom- 
passes the very general case with arbitrary number of electron 
types and arbitrary Coulomb interaction strength. Inside the 
electron bubbles summation goes over all electron branches 



In order to evaluate the decay rate and the electron 
cooling rate we need to calculate the imaginary part of 
the phonon self energy S. It is represented by the dia- 
grams shown in a Fig. S3. The second diagram is impor- 
tant when the screening is essentially dynamic. Below we 
demonstrate few particular examples how the diagram- 
matic description works. 



A. Imperfect screening 

In this Subsection we consider the case of identical 
spectral branches, but assuming now that screening is 
incomplete and electron density variations are allowed. 
Then the full e-ph vertex given by the sum of Eq.(S31) 
and Eq.(S32) is diagonal: 



l + N f vV (q) 



(S34) 



where Nf is a number of identical electron branches. 
Typically it is equal to Nf = 2N V where N v is the num- 
ber of identical valleys in a semiconductor. For example, 
N v = 6 for a bulk silicon or Nf = 4 for graphene. The 



first diagram from the Fig.S3c thus gives 

2 



Si = N f x 



1 + 2uNfV (q) 



q 2 N f v- 



-iuj + Dq 2 



(S35) 



The second diagram turns out to be crucial for a dynam- 
ical screening regime, when lo > Dq 2 : 



S 2 = 



1 + 2uN f V (q) 
1 



NfV 



—iu + Dq 2 



2 

(S36) 



V 1 {q)+N fV Dq 2 /(-iLO + Dq 2 ), 
Summing up these contributions we obtain X = Ei +£2: 



1 + NjvVoia) 



(S37) 



, lU J (lU + Dq 2 + (Nj V V (q)Dq 2 ) 



uj 2 + (Dq 2 ) 2 (l + NfvV (q)y 

The phonon decay rate is determined by the imaginary 
part of E: 

_ t T 2 NjvDq 2 

Tph ' ML ~ P^vl + iDqYil + NfvVviq)) 2 ' 
and the corresponding enhancement factor is 

(T/p F v F ) 2 



(S38) 



Tc(a) = 1 + 



1 



Cl {Vs/vp) 2 + de 2 (q 2 l 2 )(l + 2vV {q)) 2 

(S39) 

These results coincide with the ones in the main text for 
Nf = 2. Thus, we have shown equivalence of the ap- 
proach using macroscopic kinetic equation and diagram- 
matic technique. 



B. Several spectral branches 

Now we switch to the situation of complete screening of 
Coulomb interaction, so the electro-neutrality condition 
is obeyed exactly. In such a case the difference in the 
coupling constants Ti corresponding to different electron 
branches is crucial. We consider here the simplest case of 
two spectral branches. Then the full screened e-ph vertex 
given by the sum of Eqs.(S31) and (S32) is traceless: 



r s = 



ri-r 2 
2 



r,-r. 



The diagrams shown in Fig.S3c give 



2 x 



ri-r a 



2 

q v- 



(S40) 



(S41) 



— iu + Dq 2 

Thus, the phonon decay rate due to the Mandelstam- 
Lcontovich mechanism is 
1 



ph 



(rd-r 2 ) 2 



pm£(c«;,g)] 



uj—sq 



(S42) 



vDq 2 



2p n 



(Dq) 2 ' 
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and the total decay rate is enhanced (with respect to the 
classical Pippard result) by the factor 



Hq) 



1 

2q 



rx-Tg 

p F v F 



(Dqf 



(S43) 



If the asymmetry in spectral branches arise due to the 
Zeeman splitting, then T\ — T2 = (2/d e )g/iH: 



F H {q,h) = 1 



„2 u2 



did vl + (Dq) 2 



(S44) 



where h = g/iH/2£p <C 1 is dimcnsionless magnetic field 
(wc neglect here ^.-dependencies of other parameters, like 
difference D\ ^ D2, which is negligible at h <C 1). 

Note that the effect is absent in zero magnetic field, as 
the time-reversal symmetry does not allow spin density 
fluctuations to be excited by phonons. 




FIG. S4: Diffuson self energy 3, V = (l/27n/r)(l-3) . Cap- 
ital letters A and R stand for advanced and retarded electron 
Green functions respectively. 

sume a relatively weak spin-orbit (SO) interaction lead- 
ing to the spin-orbit band splitting Aso *C Ah = g/isH. 
For definiteness we consider the Rashba-typc SO coupling 
with the spin-dependent part of the Hamiltonian being 
equal to 



Hh+so = — tto x H Tr-{^xn y - a y n x ), (S47) 



A 



so 



C. Spin-orbit coupling and interbranch scattering 

In this Subsection we revisit the case of the two in- 
equivalent branches of electron spectrum and consider 
the (previously neglected) role of the interbranch scat- 
tering, using the Zeeman-splitting as an example. Qual- 
itatively, the spin- flip scattering with a rate r" 1 leads 
to non-conservation of the total spin and thus limits the 
magnitude of any effect which is related to slow spin dif- 
fusion. Formally it is described by the modification of 
the spin diffusion propagator: 



Dq 2 



Dq 2 



-ioj + Dq 2 



Dq 2 



which leads to the replacement of Eq.(S44) by 
4 v 2 F q 2 h 2 



F H {q,h) = 1 



d 2 e d u 2 + (Dq 2 +t^) 2 ' 



(S45) 



(S46) 



Below we calculate r r for the case of 2D electron gas 
with the Zeeman splitting induced by magnetic field ap- 
plied in the 2DEG plane in the x-axis direction. We as- 



where n = p/|p| is a unit vector in the direction of mo- 
mentum. The elastic scattering time turns out to be 
equal for both quasiparticlc branches(in the absence of 
electron- hole asymmetry): 



G R (e,p) E e _^ A (n)/2 + V2r' 



(S48) 



P ± = - \1 ± — — v v — ) (S49) 



A(n) 



A 2 



A 2 -2A H A so n x . 



(S50) 



where G R is the retarded electron Green's function. In 
order to find the relaxation rate, we evaluate the diffu- 
son self energy for zero frequency and momentum (lu = 
0, q = 0), Fig.S4: 



1 



dp 

2lTVT J (27r)' 



G A (0,p)®G R (0,p) 



(S51) 



A simple calculation in a manner similar to that of Rcf.7 
leads to the following result for the diffuson self energy 
at q — lu = 0: 



irA 



H 



H 



2A% 



l i3 



-si 



2 (3 + r^A|) 
H {l + T 2 A 2 H f 



s , , 



Air 3 A 3 H 
(l + r 2 A 2 H y 



r 



(S52) 



with S = (I®?— <r<g>l)/2 being the total spin of electron- 
hole pair. We are interested in the S x = subspace only 
as it hosts two eigenvalues of our interest. Naturally, the 
singlet mode(5 = 0) corresponding to the charge den- 
sity propagation remains unaffected, Ss = o = 1, while the 



triplet mode(S* = 1,S X = 0) representing a spin density 
diffusion does decay: 



A 2 
^so 



9 A 2 



(S53) 
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leading to the following result for the spin decay rate 



1 



T7 1 =T^ = - 



A 2 
^so 



2A% 



rr- 1 « r- 1 



(S54) 



In the course of derivation of Eq.(S53) we used an iden- 
tity (S = 1,S X = 0\S 2 \S = 1,S X = 0) = 1. We empha- 
size that Eq.(S54) was derived for weak SO interaction, 



IV. ANGULAR DEPENDENCE OF ULTRASONIC 
ATTENUATION 



We start here by the quasi-2D case when the thick- 
ness of a semiconductor film is much larger than the 
Fermi wavelength but still smaller than the phonon wave- 
length, Xp <C b <C X p h- In this case electron diffusion is 
two-dimensional and only the component of phonon mo- 
mentum parallel to the plane enters into the diffusion 



propagator where a replacement Dq 2 /(- 



Dq 2 ) 



Dqi) should be made. The result is that 



Eqs.(S44,S39) should be replaced by 



■Fc(«) = 1 



cf 1 {T/p F v F ) 2 sw 2 9 



{v s /v F ) 2 
2 



i-h 2 sin 



(2/d e ) 2 g 2 D (e 2 /ehv F ) 
2 a 



$55) 



Q v 2 + {Dq) 2 sin 4 6 



(S56) 



where the last equation is given for relatively strong 
screening vV (q) ^> 1. 

The true 2D case, however, should be discussed spe- 
cially. While the result for the case of imperfect screening 
at T S> p F v F is identical to Eq.(S55), the magnetic-field 
induced effect (arising from the momentum-dependent 
part of the electron-phonon vertex) may behave differ- 
ently. 

For a sufficiently thin film electron motion in the direc- 
tion perpendicular to the plane is fully quantized, thus 
the expression for electron-phonon vertex becomes 



1 



1,2 = r 6s divuH (p z ) z (iq z u z ) 

m 



(S57) 



■(P 2 z )z COS 2 



2m 



sin 2 6 



divu, 



where 



^=e F ± 9 , B H-^ 
2m 2m 



(S58) 



with ± corresponding to the spin-up and spin-down elec- 
trons. Here we have taken an average (-) z over the ground 
state corresponding to the motion in the perpendicular 



180° 




210° 



FIG. S5: Angular dependence of ultrasonic attenuation for 
magnetic-field-controlled case, Eq.S61. Red and blue curves 
represent (7y = 0, 7x = 0) and (7y = 0, 7x = 1) respectively. 
Both plots are given for InSb sample from the main text at 
frequency cj = 200 MHz 



direction and is the momentum of an in-plane mo- 
tion. We assume that the matrix element (p 2 ) z does not 
depend on the spin degree of freedom as well as on the 
electron density. Thereby we disregard any possible or- 
bital effects of magnetic field and consider the Zeeman 
interaction only which results in two different momenta 



P||,t 



and 



corresponding to the in-plane motion for 



the up and the down spin projections. The momentum- 
independent component of the vertex does not contribute 
to the magnetic-field-controlled relaxation under the con- 
dition of perfect screening (according to Eq.(S32) it is 
screened out completely). Eqs.(S43) explicitly implies 
that only asymmetric part of the vertex contributes: 



I\ — T 2 = g^H sin 2 9 divu 



(S59) 



However, in real 2D electron systems, such as het- 
erostructures, the lattice strain also affects the quasi- 
particle mass, generating the additional momentum- 
dependent contributions: 

Sm^ 1 = to _1 (7_l cos 2 9 + j\\ sin 2 9) divu, (S60) 
Ti -T 2 = gpsU x [(1 +7y) sin 2 6 + ^±_ cos 2 9] divu. 

Thus, the magnetic-field-controlled enhancement be- 
comes equal to 



T H {q, h) = 1 



(S61) 



(Dq) 2 sin 4 9 



[(1 + 7||) sin 2 6 + 7_l cos 2 6] 2 



The results (S55,S56,S61) show the angular dependence 
of ultrasonic attenuation which may exhibit a character- 
istic cross-like pattern exemplified in Fig.S5. 

There are two additional issues which should be ad- 
dressed to make the above analysis really quantitative: 
(i) in a general case the incident longitudinal(transverse) 
acoustic wave reflected off the free surface produces both 
longitudinal and transverse reflected waves, (ii) in the 
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true 2D case diffusion modes could be generated by trans- 
verse phonons as well. However, these effects do not seem 
to lead to any qualitative change of our results and we 
will postpone the corresponding studies for the future. 



V. ELECTRON-PHONON HEAT FLOW 

In this Subsection we use previously obtained results 
for the phonon decay rate, Eqs.(S55-S56), to derive an 
expression for the electron-phonon heat flow in a true 2D 
electron gas structure. We start by the spin density dif- 
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fusion effects. At the lowest temperatures T ^ Tjj 



i) 



FIG. S6: Temperature dependence of electron-phonon heat 
flow enhanced by the spin-diffusion. Plotted is the ratio of 
the spin-diffusion part Jh(T) and the local part Eq.(51). 



hv 2 /D, the enhancement does not depend on tempera- In the wide temperature region < T < where 



ture and angle, being just a numerical factor: 



J 



H 



9U 



2p 



F F 



63irh 6 p n 



(S62) 



where A x = (48(1 + -f\\) 2 + 16(1 + 7||)ti + 67J) and 
.<?□ = ^fI is the dimensionless conductance of the 2DEG. 
For higher temperatures, T 3> T^, the enhancement 
factor behaves as T oc T~ 2 . However, the resulting ex- 
pression for the e-ph heat flow depends significantly on 
the angular structure of the vertex(here A2 = (2/3) (1 + 
7||)(1 +7|| +7-0): 



H 



P F h 2 



go l5ir 3 H 6 p m v 3 s 



A2 + 7I In T/T^ ) T 4 (S63) 



If the normal strain does not alter electron mass (71 = 0) 
the angular dependence J-(q, 9) does not lead to dom- 
inance of small 9 in the corresponding integral. Then 
the heat flow is a pure power-law J oc T 4 . Otherwise, 
the shape of F{q,ff) is rather peculiar (see Fig.S5) and 
an additional In T factor appears due to the contribution 



of small angles 9 H ~ Jt^/T. Note that the T 4 InT 
behavior is slower than T 6 and at temperatures 



T > TP = 



h 2 v% 



H D 



(S64) 



the effect of spin fluctuations is smaller than the H- 
independent contribution proportional to T 6 . In the 
above Eq.(S64) we denote as Tc the i?-independent en- 
hancement factor Eq.(39) averaged over angles 9. The 
spin-orbit interaction suppresses the effect of spin fluctu- 
ations at low temperatures: 



T < T 



(so) _ 



hv. 



H VDtso ' 
so that the condition for this effect to be observed is 



T, 



(SO) 



H 



(2) 
H ■ 



The temperature dependence of the e-ph heat flow in the 



most favorable case 

T (so, < T w is ghown in Fig gg 



(i) 



~ 20 mK, T^ 2) ~ 500 mK, the outgoing heat flow 
Jh(T) oc T 4 \n(T/T^ ] ). The parameters of the 2D system 
were: electron density n = 10 11 cm -2 , electron effective mass 
m = 0.1 mo, (mo is the free electron mass), the Fermi mo- 
mentum 1/feir = 12nm, the Fermi energy Ef = 2.4meV, 
the Lande g-factor g — 5, the spin-orbit band splitting 
A_r = 0.002 Ef (corresponding to the Dresselhaus splitting 
in InP at the corresponding density); the anisotropy param- 
eters in Eq.(47) are 711 = 0, 7x = 1. The magnetic field 
H — 8T corresponds to h = 0.5. 



For the effects of charge diffusion on the e-ph heat flow 
in the 2D electron system with no additional screening 
of interactions^ = const) the situation is similar to that 
of Eq.(S62): 

j 1 tpf ( ehv P \ 2 ( r \ 2 r6; 

g a 63nh e p m v^\ e 2 / \p F VF J 

where the effective dielectric constant is an arithmetic 
mean of the dielectric constants of the media on both 
sides of the 2D system: e= (ex+£2)/2. The enhancement 
is thus reduced to a temperature-independent factor. 

The result is the same for the geometry with ad- 
ditional screening by a metallic gate at temperatures 
T 2> hv s /b with b being the distance between the 2D 
electron plane and the gate. At lower temperature T 
hv s /b the effective dielectric permittivity is q-dependent: 
e(q) = e(cothqb + l)/2 w e/2qb which transforms 
Coulomb interaction into a short-range one, Vo(q) = 
4ire 2 b. The behavior becomes similar to Eqs.(S63) how- 
ever, with the crossover temperature replaced by 
T c 1] - (l/g n k F b){v s /v F ) 2 {ehv F /e 2 )E F : 



Jc 2 



PF 



g a 4807r 3 S 6 pmwf 
2 / n \ 2 



(S66) 



kp 



(S67) 



Jc 2 D + gate — 9U 



r 1 



T< T, 



(i) 



c ' 



m v s 

Similarly to (S63), Eq.(S66) contains InT coming from 
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the angles 9 C ~ jTg'/T. 

We also note that Eq.(S66),(S67) represent only the 
Mandclstam-Leontovich contribution due to charge dif- 
fusion. At high enough temperatures 



T > T, 



(2) 



c 



T (i) (Vp 



this contribution is smaller than the J^/ arising from 
local processes corresponding to the transverse phonon 
Pippard's ultrasound attenuation, Eq.(l) in the paper: 



■/<) 



9a 



Vp 



1267rS, 6 p„ 



(S68) 



Thus for temperatures T > the T 6 law is restored. 

We should note that in real heterostructures the de- 
formation potential T is in general anisotropic. However, 
this fact does not lead to any profound changes like sup- 
pression of logarithmic behavior oc InT in Eq.(S66). This 
would require a highly anisotropic deformation potential 
r oc sin 2 9 which we do not expect. 



VI. A PARTICULAR EXAMPLE: THIN FILM OF 

InSb. 

Here we consider enhancement on ultrasound atten- 
uation in an InSb thin film with electron density n = 
10 11 cm -2 , and thickness d > 10 — 20nm on the SiC-2 
substrate. We will consider the spin effect of parallel 
magnetic field applied to the film, so its thickness d is 
chosen to be relatively small (slightly larger than Fermi 
wavelength) in order to avoid orbital effects of magnetic 
field. Phonon wavevector q is parallel to the 2DEG plane. 
Effective mass of InSb is equal to m = 0.014mo [8], 
spin-orbital band splitting is A50 ~ O.llmeV [9, 10] 
and the electron mean-free-path is supposed to be rel- 
atively long, ppl = 50. At such parameters the Fermi 
energy ep = 0.02 eV, while the deformation potential 
r = 13.12 cV[ll]. A SiC>2 substrate is characterized by 



V s = 6 X 10 5 cm/s and the effective dielectric constant 
e = (l+3.9)/2 « 2.5. Finally, we use Eqs.(S46) and (S54) 
with the specified sample/material parameters. The re- 
sulting plots are given in the Fig.l of the main text. 



VII. A GENERAL EXPRESSION FOR PHONON 
DECAY RATE: A PHENOMENOLOGICAL 
DERIVATION 

Here we derive a general expression for ultrasonic at- 
tenuation using a phenomenological approach of diffusive 
electron transport. We start by the system of equations 



' <9tn (i) + divj (i) =0, 



J V {r - r') Sn u) {r') + T w divu 



(S69) 



Fourier transforming the set we get 

(-»w + D^q 2 )n^ = - K ®4>U® 
0-M = V (q) 5^n^(q) + (tq • u)r« 



(S70) 



Due to Coulomb interaction the solution for i-th branch 
depends on the dynamics of total density. Thus, the so- 
lution is 



(S71) 



where IL^'(u>,q) = q 2 / (—iuj + q 2 ) is a response 
function, = (iq • u)r (i) , 



(S72) 



Here <I> C describes dynamic Coulomb counteraction. To 
obtain the phonon decay rate we have to evaluate the 
dissipation power in a unit volume 



ft = |ERe (j«) -P*W) = IJ2 Re - * c (w, ? ))II«($W - $o(w,?))*] 

i i 

= I Elm [($«-$ c ( w , g ))n(«)($W-$ c ( w , ? ))*] . 

i 

1 



(S73) 



and the energy of acoustic wave: 



E w = — to u m 



Finally, for the phonon decay rate we get 



r", 1 = ^ = ^E Im fr«nW(rWrl (S75) 



(S74) 



PrnU 



^rWim [n^l (rW)* (S76) 
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where T s = ($M — q))/ divu can be considered as 

a dynamically screened vertex. 

A. Multiple branches 



Coulomb counteraction term takes the form 



J2 nw$w 



(S77) 



We analyze here the case of two quasiparticle branches 
and very strong bare Coulomb potential, V(q) — > oo. The 



that exactly fixes total density <5(ni + n 2 ) = 0. Thus the 
phonon decay rate becomes equal to 



ph 



-Im 



(Tt-Ta)^ (Ij - r 2 )n* 

ni + n 2 1 (ni + n 2 )* 



(ra-ronx (Ta-ron? 
-n 2 



ni + n 2 



(it + it,)* 



(ri - r 2 



-im [n^ 1 + n 2 



(S78) 



Using also the fact that k = vD we obtain: 



ph 



(]?! - T 2 ) 2 is*D*q 2 
Pm v 2 + {D*qy 



(S79) 



where v s = uj/q is the sound velocity and = {y± 1 + 
^a" 1 )" 1 , £>* = v* l {{viD 1 )- 1 + {v 2 D 2 )- 1 )- 1 are the effec- 
tive density of states and diffusion coefficient respectively. 
To obtain the total ultrasonic attenuation we also have 
to take into account the PIC result: 



(o) 



Cj 



Pf 

2p m 



'y\D\ + v 2 D 2 )q 



(S80) 



This equation coincides with Eq.(l) of the main text if 
both electron branches are identical V\ = v 2 , -Di = D 2 . 
Thus, for the total attenuation rate we obtain: 



r£=c i ^(viD 1 + v a D 2 )<t > 1 (Tl ~ 



v*D*q 2 



2pn 



and the enhancement factor T = T pl ^ / (t^) 1 -^ is 



{D,q) 2 
(S81) 



ci V £f 



viD x 



v 2 D 2 V 2 S + (D*q) 2 
(S82) 

An important particular case is that of the Zccman split- 
ting by an in-plane magnetic field H for a 2D electron 
system (d e = 2). Here an asymmetry appears between 
the spin-up and spin-down electron branches (while v\ = 
vi = v): e m = e F ± pH/2, D m) = D{1 ± i iH/2e F ). 
However, the most important asymmetry is the one in 
the vertex T that arises from the momentum-dependent 
part of the electron-phonon coupling: 



(dT/de F )nH 



(S83) 



In the simplest model, where the only effect of the strain 
upon electron spectrum is its overall shift, the density- 
dependent contribution arises from the stress ip a vp) F s 



of electron liquid only: 

T(p F ) = T +p F v F /2 -> dT/de F = 1 (S84) 
Introducing the dimcnsionless magnetic field h = 
p,H/2e F , we arrive at the result 



T H (q,h) = 1 + 8 



v 2 F h 2 (l-h 2 ) 
v 2 + (Dq(l - h 2 )f 



(S85) 



Finally we discus the effect of inter-branch scattering. In 
fact, it modifies the response function 



D,q 2 



D,q 2 



-iuj + D*q 2 —iui + D*q 2 + t 3 , 



(S86) 



where tso is the characteristic time of inter-branch mix- 
ing so labeled in analogy with the spin-orbit mixing of 
electron branches with different spin projections. Tt lim- 
its the diffusion enhancement factor 



D«q 2 



D*q 2 



L0 2 + (D lt q 2 ) 2 L0 2 + {D,q 2 + l/T ao ) 2 
and the final result becomes equal to 

q 2 v 2 h 2 (l-h 2 ) 



F H {q,h) = 1 + 8 



(qv s ) 2 + (Dq 2 (l-h 2 ) + l/r so ) 2 



(S87) 



(S88) 



B. Imperfect screening 



Another case of interest is the case of two quasiparti- 



cle branches with identical parameters 1\( 2 



(2) 



^1(2) 



v i -Di(2) = D but finite strength of Coulomb interaction. 
In this case no asymmetry is present and thus asym- 
metric electron modes cannot be excited. However, fi- 
nite Coulomb interaction and incomplete screening al- 
lows density fluctuations which diffusive relaxation leads 
to the enhancement of ultrasound attenuation and the 
e-ph energy flow: 
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(S89) 



(r/p F v F ) 2 



q d2(w 5 /M 2 + (3 2 « 2 )(l + 2^Vb(g)) s 



(S90) 



For a 2D geometry and Coulomb interaction Vo(q) = 
2ire 2 /eq which is still relatively strong, vV 3> 1, 
vV(q)ql 3* s / w f, the result acquires a simple, frequency- 
independent form: 



Fc 2r> = 1 



1 / ehvF\ 
ci9 2 n \ e 2 J 



\p F V F J 



(S91) 



where g\j is dimensionless conductance per square in e /h 
units. The expression is valid also for a quasi-2D sample, 
when the phonon wavelength A is much larger than the 
width of quasi-2D system d s . 

In the most general case, the bare Coulomb potential 
Vq (q) acting between conduction electrons can be written 
in terms of some dispersive dielectric response s(q), as 
Vo(q) = 2-Ke 2 /qe(q). Therefore Eq.(S91) can be used 
in order to extract s(q) dependence from the measured 
phonon relaxation rate. 

An important special case is presented by a 2D elec- 
tron gas with a metal gate placed nearby which addition- 
ally screens electron-electron interaction. Here e(q) = 
e/(l - e~ 2qb ) and V (q) = {2ne 2 /eq){l - er 2c > b ), where b 
is the distance between electron plane and gate parallel 
to it. As long as phonon wavelength is shorter than dis- 
tance b, the presence of the gate may be ignored, while 
for short wavelengths, when qb -C 1, Coulomb interac- 
tion becomes effectively-short-range, Vo(q) = 47re 2 6/e. 
Exactly for this region of low frequencies (temperatures) 
enhancement factor T is 



C2L)+gate 



1 



16(T/ PF v F y< 



i(v 8 /v F ) 2 + (g 2 Z 2 )(87we 2 &/e) 2 



(S92) 



We see that the crossover between T 
behavior emerges at 



h J_J_ 



ga k F b \v F J \ e 



and T 



E F . 



This equation can be used for an experimental determi- 
nation of the background dielectric constant e. 
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